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Abstract 

In this paper we study symmetry reductions of a class of nonlinear 
third order partial differential equations 

Ut — eUxxt + 2k.Ux = UUxxx + aUUx + Pu^Uxx , (1) 

where e, k, a and (i are arbitrary constants. Three special cases of equa- 
tion (1) have appeared in the literature, up to some rescalings. In each 
case the equation has admitted unusual travelling wave solutions: the 
Fornberg-Whitham equation, for the parameters e = 1, a = —1, /3 = 3 
and K = i, admits a wave of greatest height, as a peaked limiting form 
of the travelling wave solution; the Rosenau-Hyman equation, for the pa- 
rameters e = 0, a! = l, /3 = 3 and k = 0, admits a "compacton" solitary 
wave solution; and the Fuchssteiner-Fokas-Camassa-Holm equation, for 
the parameters e = 1, a = —3 and (3 — 2, has a "peakon" solitary wave 
solution. 

A catalogue of symmetry reductions for equation (1) is obtained using 
the classical Lie method and the nonclassical method due to Bluman and 
Cole. 



1 Introduction 

In this paper we are concerned with symmetry reductions of the nonlinear third 
order partial differential equation given by 

A = Uf — eUxxt + 2KWa; — UUxxx — OLUUx — (3UxUxx = 0, (1-1) 

where e, k, a and (3 are arbitrary constants. Three special cases of {^) have 
appeared recently in the literature. Up to some rescalings, these are: (i), the 
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Fornberg-Whitham equation [28,58,59], for the parameters e — 1, a ^ — 1, 
P — 3 and k — ^, (ii), the Rosenau-Hyman equation [54] for the parameters 
e — Q,a — 1, f3 — 3 and k — 0, and (iii), the Fuchssteiner-Fokas-Camassa-Holm 
equation [9,10,25,27] for the parameters e = 1, a = — 1 and /3 = 2. 
The Fornberg-Whitham (FW) equation 

Ut — Uxxt + Ux = UUxxx — UUx + 3u (1.2) 

was used to look at quahtative behaviours of wave-breaking [58] . It admits a 
wave of greatest height, as a peaked limiting form of the travelling wave solution 
[28], 

u{x,t) — Aexp (— ^[a; — |t|) , 

where A is an arbitrary constant. 

The Rosenau-Hyman (RH) equation 

Ut = UUxxx + UUx + 3UxUxx- (1-3) 

models the effect of nonlinear dispersion in the formation of patterns in hquid 
drops [54]. It also has an unusual solitary wave solution, known as a "com- 
pacton" , 

-fccos2{i(a;-ci)}, if [a; - ct| < 27r, 
0, if \x - ct\ > 2ti. 



u{x, t) 



These waves interact producing a ripple of low amplitude compacton- 
anticompacton pairs. 

The Fuchssteiner-Fokas-Camassa-Holm (FFCH) equation 

Ut - Uxxt + "^KUx = UUxxx - "iuUx + 2UxUxx, (1-4) 

first arose in the work of Fuchssteiner and Fokas [25,27] using a bi-Hamiltonian 
approach; we remark that it is only implicitly written in [27] — see equations 
(26e) and (30) in this paper — though is explicitly written down in [25]. It has 
recently been rederived by Camassa and Holm [9] from physical considerations 
as a model for dispersive shallow water waves. In the case k = 0, it admits an 
unusual solitary wave solution 

u{x^ t) — Aexp {—\x — ct\) , 

where A and c are arbitrary constants, which is called a "peakon" . A Lax-pair 
[9] and bi-Hamiltonian structure [27] have been found for the FFCH equation 



(1.4) and so it appears to be completely integrable. Recently the FFCH equation 
(1.4) has attracted considerable attention. In addition to the aforementioned, 
other studies include [10,21,22,23,24,26,32,41,46]. 

The FFCH equation (|l.4|) may be thought of as an integrable modification 



of the regularized long wave (RLW) equation [7,47] 

Uxxt + UUx - Ut- Ux ^ 0, (1-5) 
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sometimes known as the Benjamin-Bona-Mahoney equation. However, in con- 
trast to (|1.4D, the RLW equation (1.5) is thought not to be solvable by inverse 
scattering (cf., [42]); its solitary wave solutions interact inelastically (cf., [37]) 
and only has finitely many local conservation laws [45]. However physically 
it has more desirable properties than the celebrated Korteweg-de Vries (KdV) 
equation 

Ut + U:cxx + GuUx = 0, (1.6) 

which was the first equation to be solved by inverse scattering [31]. We remark 
that two other integrable variants of the RLW equation (1.5) are 



Uxxt + 2uut - Uxd„ ^Ut - Ut - Ux ^ 0, 



(1.7) 



where (9^ ^/) {x) = f{y) dy, which was introduced by Ablowitz, Kaup, 
Newell and Segur [2], and 



Uxxt + UUt 



Uxd^ ut 



Ut~Ux ^ 0, 



(1.8) 



which was discussed by Hirota and Satsuma [34]. We also note that (|L 
with K = i, (|r5|), (|r^) and (^) all have the same linear dispersion relation 
(jj{k) = —k/(l + fc^) for the complex exponential uix, t) ^ exp{i[fca; + a;(fc)i]}. 

Recently, Gilson and Pickering [32] have shown that no equation in the 
entire class of equations (1.1) will satisfy the necessary conditions of either the 
Painleve PDE test due to Weiss, Tabor and Carnevale [57] or the Painleve ODE 
test due to Ablowitz, Ramani and Segur [ 3,4] t o be solvable by inverse scattering. 
However, the integrable FFCH equation (1^) doe s po ssess the "weak Painleve" 
property (cf., [49,50]), as does the FW equation (|l^ ). 

All these special travelling wave solutions are essentially exponential solu- 
tions, or sums of exponential solutions, and thus would suggest some sort of 
linearity in the differ enti al equation. This is discussed by Gilson and Pickering 
[32], who show that (1.1), with a 7^ and (3{1 + (3) ^Q, can be written as 



{I3ux + udx + edt) {uxx 



2k/ = 0, 



(1.9) 



where dx = d/dx, dt = d/dt and ^j? = —a/{l+(3), provided that ea + (3+l = 0, 
which includes the FFCH equation (1^). For the traveling wave reduction, 

u = w{z), z = X — ct, 

the resulting ordinary differential equation is 

(Zk — c)w + ecw — WW — aww — pw w — 

where ' = d/dz, which also may be factorised as 



(1.10) 



Pw' + {w 



{w" - fi'^w + 7) = 0, 



(1.11) 



provided that 



1 + /3' 



^(1 + l3)-i - 2k(1 +l3) + c{l+l3 + ae) = 0. 
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This includes all three special cases (|1.2|)-(1.4); since is strictly non-zero 



in these three cases then a suitable 7 can always be found. 



Furthermore, ii 1 + P + ae — and e 7^ 0, then (1.1) with n — possesses 
the "peakon" solution 

u{x,t) = Aexp (^-e"^/^|a; - ct\j , 
where A and c are arbitrary constants. More generally, if a/{l + /?) < 0, 



1 + /3 + ae ^ and k ^ 0, then (1.1) possesses the solution 



. N . if a V^^ , ,1 2(1 + /3)k 

.(x,0=Aexp|-(^-^j \x-ct\j, c=^^^^, 

where A is an arbitrary constant. If a/(l + /3) > 0, /3 —1 and a/3 7^ 0, then 
( p~l| ) possesses the "compacton" solution 

^ 2|2(l + ffl.-a + g.a.)c| j , ,^ _ 

where c is an arbitrary constant. 

The classical method for finding symmetry reductions of partial differential 
equations is the Lie group method of infinitesimal transformations. As this 
method is entirely algorithmic, though often both tedious and virtually unman- 
ageable manually, symbolic manipulation programs have been developed to aid 
the calculations. An excellent survey of the different packages available and a 
description of their strengths and applications is given by Hereman [33] (see also 
his contribution in this volume). In this paper we use the MACSYMA package 
symmgrp . max [11] to calculate the determining equations. 

In recent years the nonclassical method due to Bluman and Cole [8] (in the 
sequel referred to as the "nonclassical method"), sometimes referred to as the 
"method of partial symmetries of the first type" [56], or the "method of con- 
ditional symmetries" [35] , and the direct method due to Clarkson and Kruskal 
[15] have been used to generate many new symmetry reductions and exact solu- 
tions for several physically significant partial differential equations that arc not 
obtainable using the classical Lie method (cf., [13] and the references therein). 
The nonclassical method is a generalization of the classical Lie method, whereas 
the direct method is an ansatz-based approach which involves no group theo- 
retic techniques. Nucci and Clarkson [43] showed that for the Fitzhugh-Nagumo 
equation the nonclassical method is more general than the direct method, since 
they demonstrated the existence of a solution of the Fitzhugh-Nagumo equa- 
tion, obtainable using the nonclassical method but not using the direct method. 
Subsequently Giver [45] (see also [6,48]) has proved the general result that for 
a scalar equation, every reduction obtainable using the direct method is also 
obtainable using the nonclassical method. Consequently we use the nonclassical 
method in this paper rather than the direct method. 
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Symmetry reductions and exact solutions have several different important 
applications in the context of differential equations. Since solutions of par- 
tial differential equations asymptotically tend to solutions of lower-dimensional 
equations obtained by symmetry reduction, some of these special solutions will 
illustrate important physical phenomena. In particular, exact solutions arising 
from symmetry methods can often be used effectively to study properties such as 
asymptotics and "blow-up" (cf., [29,30]). Furthermore, explicit solutions (such 
as those found by symmetry methods) can play an important role in the de- 
sign and testing of numerical integrators; these solutions provide an important 
practical check on the accuracy and reliability of such integrators (cf., [5,55]). 

Classical symmetries of differential equations are found in practice by a two- 
step process. The first involves finding the determining equations for the in- 
finitesimals of the group action. These determining equations form an overde- 
termined, linear system of partial differential equations. The second step in- 
volves integrating this system. The first step is entirely algorithmic, and has 
been implemented in all the commercial symbolic manipulation languages (cf., 
[33]). The second step involves heuristic integration procedures which have been 
implemented in some symbolic manipulation programs and are largely success- 
ful, though not infallible. Commonly, the ovcrdctermined systems to be solved 
are simple, and heuristic integration is both fast and effective. However, there 
are three areas where heuristics can break down (cf., [39] for further details and 
examples). 

1. Arbitrary param,6ters and functions. If the partial differential equation 
whose symmetries are sought involves arbitrary parameters, such as (1.1) 
or more generally, arbitrary functions, heuristics yield usually the general 
solution, and miss those special cases of the parameters and arbitrary 
functions where additional symmetries exist. 

2. Termination. Heuristic algorithms are not guaranteed to terminate, and 
may become trapped in infinite loops for some examples. 

3. Too difficult to solve. The system may not be solvable by the heuristic. 
The heuristic will then attempt to represent the general solution in terms 
of functions satisfying certain conditions, but may give up before a useful 

representation is obtained. 

These problems are addressed by use differential Grobner bases (DGBs) which 
we describe below. 

The method used to find solutions of the determining equations in the non- 
classical method is that of DGBs, defined to be a basis Q of the differential ideal 
generated by the system such that every member of the ideal pseudo-reduces 
to zero with respect to fi. This method provides a systematic framework for 
finding integrability and compatibility conditions of an overdetermined system 
of partial differential equations. It avoids the problems of infinite loops in re- 
duction processes and yields, as far as is currently possible, a "triangulation" of 
the system from which the solution set can be derived more easily [16,40,51,52]. 
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In a sense, a DGB provides the maximum amount of information possible using 
elementary differential and algebraic processes in finite time. 

In pseudo- reduction, one must, if necessary, multiply the expression being 
reduced by differential (non-constant) coefficients of the highest derivative terms 
of the reducing equation, so that the algorithms used will terminate [40]. In 
practice, such coefficients are assumed to be non-zero, and one needs to deal 
with the possibility of them being zero separately. These are called singular 
cases. 

The triangulations of the systems of determining equations for infinitesimals 
arising in the nonclassical method in this paper were all performed using the 
MAPLE package diffgrob2 [38]. This package was written specifically to han- 
dle nonlinear equations of polynomial type. All calculations are strictly 'poly- 
nomial', that is, there is no division. Implemented there are the Kolchin-Ritt 
algorithm using pseudo-reduction instead of reduction, and extra algorithms 
needed to calculate a DGB (as far as possible using the current theory), for 
those cases where the Kolchin-Ritt algorithm is not sufficient [40] . The package 
was designed to be used interactively as well as algorithmically, and much use 
is made of this fact here. It has proved useful for solving many fully nonlinear 
systems [16-19]. 

In the following sections we shall consider the cases e = and e ^ 0, when 
we set e = 1 without loss of generality, separately because the presence or lack 
of the corresponding third order term is significant. In §2 we find the classical 



Lie group of symmetries and associated reductions of (1.1). In §3 we discuss the 



nonclassical symmetries and reductions of (1.1) in the generic case. In §4 we 
consider special cases of the the nonclassical method in the so-called t = 0; in 
full generality this case generates a single equation which is considerably more 
complex than our original equation! In §5 we discuss our results. 

2 Classical symmetries 

To apply the classical method we consider the one-parameter Lie group of in- 
finitesimal transformations in (x, t, u) given by 

X* = X + e£,{x,t,u) + 0{e'^), 

t* = t + ST{x,t,u) + 0{e'^), (2.1) 
u* = u + e(p{x,t,u) + 0{e'^), 

where e is the group parameter. Then one requires that this transformation 
leaves invariant the set 

Sa = Mx,t) : A^O} (2.2) 



of solutions of (1.1). This yields an overdetermined, linear system of equations 
for the infinitesimals ^(x, t, u), t{x, t, u), 4>{x, t, u). The associated Lie algebra is 
realised by vector fields of the form 

V = £^{x, t, u)dx + t{x, t, u)dt + <j){x, t, u)du- (2.3) 
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Having determined the infinitesimals, the symmetry variables are found by solv- 
ing the characteristic equation 

da: = ^ d"" ^2 4) 

^{x,t^u) T{x,t,u) (j){x,t,u)' 

which is equivalent to solving the invariant surface condition 

■0 = ^(x, t, u)ux + t{x, t, u)ut — 4'{x, t, u) — 0. (2-5) 

The set S'a is invariant under the transformation (^]^) provided that 
pr(3)v(A)|A=o = where pr^^V is the third prolongation of the vector field 
( ^.3| ), which is given explicitly in terms of and (j) (cf. [44]). This procedure 
yields the determining equations. There are two cases to consider. 

2.1 e = 

In this case using the MACSYMA package symmgrp . max we obtain the following 
system of ten determining equations 

Tu =0, Tj; = 0, = 0, U(puuu + P(l>uu = 0, iu^cjiuu + - /3(/) = 0, 

nu - i^xU + (/) = 0, 4>xxxU + {au - 2k)4>x = 0, 
5u'^4>xxu + Pu4>xx + '2K(t> - u^ixxx + (2au^ - 4km)^^ + U^t = 0. 

(2.6) 

Next applying the reduceall algorithm in the MAPLE package dif f grob2 to 
this system yields 

(2 + /3)e.. = 0, (2 + /?)[a<,t + - 2<.t] = 0, 

£.u =0, Tx= 0, Tu = 0, 

2q;<j; + 2k^^ +£.t~ 2KTt = 0, (2 + /3)[2k(/) + [iau^ - 4ku)^j; + <(] = 0. 

This is simple enough to solve; there is no need to do the full Kochin-Ritt 
algorithm in this case. The output shows that there are three special values of 
the parameters, namely a = 0, /3 = — 2 and k = 0, and combinations thereof. 
It transpires that the special case /5 = — 2 is purely an artefact. For the three 
special cases (a) a = 0, k ^ 0, (b) a 7^ 0, k = and (c) a = k = 0, applying 
the reduceall algorithm of diffgrob2 to (2.(;) yields 

(a) a = 0, K^O ^xx^O, 6t-2<,t = 0, = 0, 

= 0, 2k^j; + it - '2KTt = 0, T„ = 0, 

2k(/) + (2au2 - 4Ku)£,a: + uit = 0. 

(b) a^O, K = 2auC, +it = 0, ^tt = 0, Cu = 0, 

Ta; = 0, Ttt = 0, T„ = 0, 

2k0 + {2au'^ - 4:Ku)ix + <t = 0. 
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(c) a = K = ^^0^ = 0, = 0, = 0, 

Tx = 0, Ttt = 0, Tu= 0, 
— 3u^a: + UTt — 0. 

Hence we obtain the following infinitesimals: 
Case 2.1(i) a 7^ and k 7^ 

^ = 2KC3t + ci, r = C3t + C2, 4'=-c^u. (2.7) 
Case 2.1(ii) a = and k 7^ 

^ = + 2k(c4 — ca)^ + ci, T = + C2, = (3c3 — 04)11. (2.8) 

Case 2.1(iii) a 7^ and k = 

C = Ci, T = C3t + C2, (f) ^ -C3U. (2.9) 

Case 2.1(iv) a = and k = 

^^cax + ci, T = C4i + C2, = (3c3 — C4)m, (2-10) 

where ci, C2, C3, C4 are arbitrary constants. 

Solving the invariant surface condition ( ^.5[ ) yields four different canonical 
reductions: 



Reduction 2.1 a and k arbitrary. If C3 = C4 = in (2.7)-(2.10) we may set 
ci — c and C2 = 1 and we obtain the travelling wave reduction 

u{x, t) — w{z), Z — X — ct, 

where w{z) satisfies 

ww'" + aww' + fiw'w" + (c - 2n)w' = 0. 

This can be integrated to yield 

ww" + l){w'f + ^aw"^ + (c - 2k)w = A, 

where A is an arbitrary constant. Multiplying this by w^^'^w' and integrating 
again yields 

+ + + (2.11) 

where B is an arbitrary constant, for (3 7^ —1,0,1. Generally if /3 7^ —1,0,1, 
then (2.11) is solvable using quadratures, though for certain special values of 
the parameters there are explicit solutions. For example (i), ii f3 — —2 or /3 = 
—3, then (2.11) is solvable in terms of Weierstrass or Jacobi elliptic functions. 
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respectively, (ii), if B = 0, then ( ^.11 ) is solvable in term of trigonometric 
functions, and (iii), if c = 2k and /? = 3, then w{z) can be expressed in terms 
of trigonometric functions via the transformation w(z) = w^/^. 
In the special cases /3 = — 1, 0, 1 we obtain the equations 

(w')^ + ctw'^ Inw + 2(2k - c)w = Bw'^ - A, 
(w')^ + "10^ + 2(2k - c)w\nw Bw - 2A, 
(w')^ + aiy^ + 2(2k - c)w = B - A\nw, 

respectively, with A and B arbitrary functions. If the coefficient of Inw in these 
equations is zero, then w{z) is expressible in terms of elementary functions, 
otherwise in terms of quadratures. 



Reduction 2.2 a ^ 0, k arbitrary. If C3 7^ in (2.7) and (2.£) we may set 
C3 = 1, ci = c and C2 — 0, without loss of generality, and obtain the reduction 



u{x, t) = w{z)t z 
where satisfies 



clnt — 2K.t, 



(2.12) 



' Pw'w" + aww' + cw' 



w = 0. 



Also if C3 = and C4 7^ in (2.8) we may set — I, ci ~ c and C2 — 0, without 
loss of generality, and obtain the same reduction ( ^.12| ). 

Reduction 2.3 a = 0, k arbitrary. If C3 7^ and C4 7^ in (2.8) and ( 2.10| ), 

we may set C3 = m + i, C4 = 1 and ci = C2 = 0, without loss of generality, and 
obtain the reduction 

u{x, t) = w{z)t^"', z = {x- 2Ki)f-™-i/3, 

where w{z) satisfies 

ww'" + fiw'w" + (m + \)zw' — 3mw — 0. 



Reduction 2.4 a — 0, k arbitrary. If C3 7^ and C4 = in (2^) and ( 2.1C ), we 

may set C3 — m, ci — 2k and C2 — 1, without loss of generality, and obtain the 
reduction 

u{x,t) = wiz)e^"^\ = (x- 2Kt)e-"*, 
where w{z) satisfies 



WW 



■ f3w'w" 



3mw = 0. 
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(2.13) 



2.2 e=l 

In this case we obtain the following system of eleven determining equations 

Tu ==0, Tx = 0, ^„ = 0, 0„„ = 0, 2(j).^u - ixx = 0, 

I3{u(l)u " + Ct) = 0) 4> + uTt - u£,x - = 0, 
3u0xij + (j^tu + - iu^xx - 2Cxt = 0, 
U4>xxu + + UTf - 3^a;W - = 0, 
U(l)xxx + 0a:xt + (aU " '2l^)4>x = 0, 

+ (201*2 - 4K)C:r + [(a + 1)m - 2K]^t = 0. 

As in the previous case, we apply the reduceall algorithm in the MAPLE 
package dif f grob2, to this system, which yields 

=0, (a + l)iu = 0, C« = 0, 

Tx =0, (a + l)Ttt =0, T„ = 0, 

2k(/« = [2k- (« + l)u]^t. 

This shows that there are two special values of the parameters, namely a = — 1 
and K = 0. For the three special cases (a) a = —1, k 7^ 0, (b) a 7^ —1, k = 
and (c) a — —1, k = 0, applying the reduceall algorithm of diffgrob2 to 
(pl3|) yields 



(a) a = — 1, K 7^ 



= 0, 6t = 0, = 0, 

Tx =0, Tt = 0, T„ = 0, 



(b) a^-1, K = ex=0, 6 = 0, C« = 0, 



0, 



1, K = 



0, 



Tx = 0, Ttt = 0, r, 

= 0, itt = 0, 

Tx = 0, Tu = 0, r„ = 
(f)- u£,t+ uTt = 0. 

Hence we obtain the following infinitesimals: 



Case 2.2(i) a 7^ -1, k 7^ 

(1 + Q:)c3t 

4=C3t + Ci, T= ^- hC2, = C3 

Case 2.2(ii) a = -1, k 7^ 

^ = C3i + Ci, T = C2, = 



(1 + a)u 
2^. 



C3- 



(2.14) 
(2.15) 
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Case 2.2(iii) -1, k = 

C = Ci, T^C3t + C2, ^ = -C3M. (2-16) 

Case 2.2(iv) a = -1, k = 

^ — Cst + Ci, T — C4t + C2, = C3 — C4W, (2-17) 

where Ci, C2, C3, C4 are arbitrary constants. 
There are four canonical reductions. 



Reduction 2.5 a and k arbitrary. If in (2.14)-(2.17) — 04, — 0, we may set 

ci = c and C2 — I without loss of generality. Thus we obtain the reduction 

u{x, t) — w{z) + c, z — X — ct, 

where w{z) satisfies 

ww'" + I3w'w" + avow' — [2k — (1 + a)c\w' . 

This can be integrated to yield 

ww" + i(/3 + l){w'f + \aw'^ = [2k - (1 + a)c]w + A, 

where A is an arbitrary constant. Then multiplying through by w^^'^w' and 
integrating again yields 

[w'f + 1^ = + + 5^1-/., (2.18) 

p + 1 p P — 1 

provided that (3 ^ —1,0,-1. Generally if j3 7^ —1,0,1, then ( ^.1S| ) is solvable 
using quadratures, though for certain special values of the parameters, there 
are explicit solutions. For example (i), if f3 = —2 or (3 = —3, then ( 2.18 ) is 
solvable in ter ms o f Weierstrass or Jacobi elliptic functions, respectively, (ii) if 
B = 0, then ( 2.18 ) is solvable in term of trigonometric functions, and (iii) if 
(1 + a)c — 2k, and /3 = 3, then w{z) can be expressed in terms of trigonometric 
functions via the transformation w{z) ~ w^/^. 

In the special cases /? = — 1, 0, 1 we obtain the following equations, 

[w'Y + 2aw'^ \nw = -2[2k - (1 + a)c^w - A + Bw'^, 
(w')^ + 2aw^ = -2[2k- (1 + a)c]wlnw - 2A + Bw"^, 
{w')^ + 2aw'^ = -2[2k - (1 + a)c\w + 2Aliiw + Bw"^, 

respectively, where A and B are arbitrary constants. If the coefficient of Inw 
in these equations is zero, then w{z) is expressible in terms of elementary 
functions, otherwise in terms of quadratures. 
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Reduction 2.6 a ^ —1, k arbitrary. If C3 7^ in (2.14), we may set C3 = 1, 

C2 = and ci = 2/tc/(l + a), without loss of generality. Thus we obtain the 
reduction 

, , wiz) + c 2k, 2Kt 
u{x,t) = ^ + ^—, z = x-— clnt, 2.19 

t 1 + a 1 + Q' 

where w{z) satisfies 

ww'" + (3w'w" - w" + avow' + (a + l)cw' + w + c = 0. (2.20) 



If C3 7^ in ( 2.16| ) we may set C3 = 1, ci = c and C2 = to obtain the reduction 
( ^.19D with K = 0. 



Reduction 2.7 a = —1, k 7^ 0. If C3 7^ in (2.15) then we set C3 = to, ci = 
and C2 = 1, without loss of generality. Thus we obtain the reduction 

u{x^t) = w{z) + mt, z = X — \mt^ , (2-21) 

where w{z) satisfies 

III , fD I II I n I n 

WW + pw w — WW — 2kw — to = U, 
which may be integrated to yield 

ww" + i(/3 - l)(w')^ - - 2kw ~mz^ a, (2.22) 
where A is an arbitrary constant. 



Reduction 2.8 a — —1, k = 0. If C3 7^ and C4 7^ in (2.17) we may set 
C3 = TO, C4 = 1, ci = c and C2 = 0, without loss of generality. Thus we obtain 
the reduction 

u{x,t) — h TO, z = X — mt — chit, (2.23) 

where w{z) satisfies 

ww'" + Pw'w" - w" - ww' + w + c^O. (2.24) 



3 Nonclassical symmetries (r ^ 0) 

In the nonclassical method one requires only the subset of given by 

^A,^ = {u{x, t) : A(m) = 0, i!{u) = 0}, (3.1) 

where S'a is defined in ( |2.2| ) and — \s the invariant surface condition ( |2.5D , 
to be invariant under the transformation ( ^.l[ ). The usual method of apply- 
ing the nonclassical method (e.g. as described in [35]), involves applying the 
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prolongation pr'^^W to the system composed of (1.1) and the invariant surface 
condition (2.5) and requiring that the resulting expressions vanish for u € S'a,?/;, 
i.e. 



pr 



(3) 



v(A)|a=o,^=o = 0, pr(^W(-i/')|A=o,i>=o = 0. 



(3.2) 

It can well known that the latter vanishes identically when -0 = without 
imposing any conditions upon ^, r and 0. To apply the method in practice 
we advocate the algorithm described in [18] for calculating the determining 
equations, which avoids difficulties arising from using differential consequences 
of the invariant surface condition (2.5). 

In the canonical case when r ^ w e se t r = 1 without loss of generality. We 
proceed by eliminating ut and Uxxt in (1.1) using the invariant surface condition 
{2.^) which yields 



^^Uxxx ^'^xxx '^^^u^x'^xx j^^x'^xx ^^u^xx ~t~ 26^2;'^^ 

+e^ti«u^ - eftynuul + 2eS,xuul - auux - 2ecj)xuUx + 2kUx 

+£^xxUx - e4>xx + - £.Ux = 0. 



(3.3) 



We note that this equation now involves the infinitesimals ^ and (p that are to 



be determined. Then we apply the classical Lie algorithm to (3.3) using the 
third prolongation pr'^^-'v and eliminating Uxxx using (3.3). It should be noted 
that the co effici ent of Uxxx is (C — ^u). Therefore, if this is zero the removal of 
Uxxx using (3.3) is invalid and so the next highest derivative term, Uxx, should 
be used instead. We note again that this has a coefficient, /3 — 3, and so that 
in the case ^ — u one needs to calculate the determining equations for the cases 
/3 ^ 3 and f3 — 3 separately. Continuing in this fashion, there is a cascade of 
cases to be considered. In the remainder of this section, we consider these cases 
in turn. First, however, we discuss the case given by e = 0. 



3.1 e = 

The first determining equation gives — 0, and substituting this into the other 
seven determining equations yields 

0«ti«" + P4>uu = 0, 3(j)xuuU + 2P(j)xu - Kxx = 0, 
34>uuu'^ + P4>uU - /30 = 0, 3(l)xuU - 3^xxU + P4>x = 0, 

<PtU - (pxxxU^ - acj)xV? + 2K(j)xU + i£,x4'u - 0^ = 0, (3.4) 

'i4>xxuU'^ - ixxxU^ + 2aixU'^ + P(t>xxU - ^£,xK.U + ?>S,S,xU 

+£_tu + 2n(t) -^0 = 0. 

It is quite straightforward to solve these equations and so we obtain the following 
infinitesimals: (a), if a ^ 



(i) i = 2K 



Cl 



t + C2 ' t + C2 ' 

(ii) e = ci, = 0, 
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and (b), if a = 

(ci + l)x + 2k(2ci - l)t + C2 



(i) C = 



5{cit + C3) ' Cli + cs' 



r-\ c x + 4:Kt + ci 

These are all equivalent to classical infinitesimals. Hence in this case there are 
no new nonclassical symmetries. 

3.2 e = 1 

As discussed in the preamble to this section, we must consider, in addition to 
the general case of the determining equations, each of the singular cases of the 
determining equations. 

Case 3.2.1 ^ ^ u. We can remove factors of — u) from the determining equa- 
tions, and we have then that = 0. Reducing the remaining eight determining 
equations with respect to this, only the last six are non-zero: 

3(t>uuu'^ - &^<i>uuu + I34>uu + 3^V«u - /3^<^« -I3<i> + P^^c + P^t = 0, 

^x(t>uU - Pi4>x + <i>4>uuU + l34>xU - ^(j)4>uu - 5^0x«W + 4^CxxM + 4>tuU 
-(t>4>u + ^t(t>u - ^4>tu - ^'^^xx + 34>xuU'^ - HxxU"^ - 2CxW - ^^xtU 

+2^Vx« + 2^x<^ - 26^x + 2Cext = 0, 

2^K(f)x - ^tU + Oi(j)xU^ - 2k^xU - a^lpxU + '2^xu^xU + ^(pxxuU - 2^(pxu4'x 

-£,xx4>xU - S,4>(t)xxu - SCx^W + '2£,£,x4> + ^^xx<l>x + £.x4'xxU - £,4>xxxU 
-^t(t> + 4''^ - (1^4'xx + (t'xxxU^ + (t>xxtU - £,(l)xxt + it<l>xx + i,(t>t = 0, 



2f3(j)xuU - ^x4>uuU - (3^xxU + 24>u(j>uuU + K^xx " H(t>xuuU ^y^uuu 
+4'4'uuuU + <j)tuuU - (t)4>uu + itft'uu - i(t>tuu + ^<t>xuuU'^ + 2£.^4'xuu 

-2^(j)u4'uu + 2£,S,x4>uu - 2P^(j)xu = 0, 

4^j.KM - 2(t>uu(l>xU - (i4>xxU - i^xx4>u - 2«;0 - 2(t>4)xuuU - 24>xtuU + C0 
+ixxtU + S^xxxyr + S^tixx - ^<l>xxuU + 2#x« - " £,£,xxt + 2^ ^x ~ C 4'xxu 

^xx^u'^ 2^CxK - ^(.xxxU - H^xU + '2'^<f>uu4x 
+^ixixx + 4:C(f>xxuU + (3^(j>xx + 2^(f)(f>xuu + 2^0„^!>x« " 2^£,x<t>xu " Ol^xU^ 
+Ol(.(.xU + 2^(pxtu + 26k - ^tU - ^xx<P - '^<t>u(l>xuU = 0. 

Reducing the fifth of these equations with respect to the fourth yields 

{13-3) [{u-O<l>u-'P + ^^x + ^t]=0. 
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If /? = 3, then one easily finds via another route that the expression in the 
second bracket is necessarily zero. The equation for 4> can be solved to give 



When this is substituted into the remaining equations we can then take 
coefficients of powers of u to be zero, and our problem is then easily solved. As 
in the e = case discussed in §3.1 above, it is quite straightforward to solve 
the resulting equations. The complete solution set is 

(a) , if a 7^ —1 

(i) ^ = ci, q> = 0, 

_ 2k Cl 2k - (1 + a)u (3.5) 

^ ^ (1 + a) < + C2 ' " (1 + a){t + C2) ' 

(b) , if a = -1 

^ = cii + C2, = ci, (3.6) 

(c) , if a = — 1 and k = 

^ = (3.7) 

(d) , if /? = -1 and a = 

= cix — 2ciKt + C2, (j) — 3ciu — 2cfx + iciKt — 2ciC2 — 2ciK, /3 7^ 0. (3.8) 



The infinitesimals (|3.5| )~ (|3.7D give rise to classical reductions, but ( |3.8| ) gives 
the following new nonclassical reduction. 



Reduction 3.1 If in ( |3.8| ), we set ci 7^ and C2 = 0, without loss of generality, 
then we obtain 

u{x,t) = w{z) exp (3cit) + C12: exp (cii) + 2k, z — {x — 2Kt — 2k/ci) exp (— Cit) , 
where w(z) satisfies 



WW 



III 



w'w" + cizw' — iciw = 0. 



Case 3.2.2 ^ = u, P ^ 3, (3 ^ 1. We generate five determining equations, the 
first of which is (j)uu = 0. Thus (/) is a linear function of u, and substituting this 
into the remaining four determining equations, we take coefficients of powers 
of u to be zero. These equations are easily solved to give = provided that 
Av = and a = —1. The invariant surface condition and ( |l.l| ) are then solved to 
give the simple exact solution 

X + Cl 

u(x,t) = — , 

t + C2 
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where Ci and C2 are arbitrary constants. 

Case 3.2.3 ^ = u, = 1. We consider here the case (?i„„ ^ 0, since taking 
'Puu = yields the same sohition as in Case 3.2.2 above. In this instance the 
remaining four determining equations are 

12k - 24>xuuU - 6au - 6m - 2^(/)„„„ - 2>(j)u4'uu - ^(l^xu ~ 2(^t„„ = 0, 

4>xu4>xxU - (t>u4>xxxU - (yAAxU - (f>(j)xuU - '2(j)x4>xx + 4'4>uu4>xx 

+ <t>tu<Pxx - '2(pu4>xu4>x + '2K(Pu4'x + 2000; " " (jxt^uf/^XXU 

-4>xxt(t^u + 4>t(t>u - 4'4>tu = 0, 

4>u4'xuuU + 4a0„U - (j)(f>l^ ~ (t>tu4>uu + 4>4>u4'uuu +&4>- Hxx 
+0t«u0« - 4:(j)xxuU + 4:(j)uU - 24>uu(l>x - H4>xuu + 2(/)^0„„ - 8k(/)„ 

-2a(j) - <i)xu4>uuU - ^(t>xtu = 0, 

(t>(l>uuU + a(j)xuU'^ - 2(i)u4>uu4>x - 24>(j)u<i)xuu - 2(j)u4>xxuU + a(f)tuU 
+2(j)(t)xu^uu - '2,K(j}(j)uu + 2(j)tu(l>xu + OicfxpuuU + 2(j)xxxU - 2(j)xU 
+2(j)1yU - 'i(l)u4>xx + 24>4>xxu - 2(l)i(l)xu - 2(pxtu4>u + H(l>u - 2K(t>tu 
-a(f)(pu - 2k0x«W + ^tuU + (t^xuU"^ + 2^xxt - 2(j)t = 0. 

Using the procedures in the package diffgrob2 with an ordering designed 
to eliminate first derivatives with respect to t, then derivatives with respect 
to X, one can obtain several equations for derivatives of </> with respect to 
u only. One can then continue to produce lower order and lower degree 
equations in the M-derivatives of 0, using repeated cross-differentiation and 
reductions. For example, the "Direct Search" procedure in the diffgrob2 
manual, [38] may be used. This process suffers from expression swell. No 
termination of this process was observed by us within the computer memory 
available, and the expressions obtained contained thousands of summands! 
One of three results appear likely. Firstly, the process terminates with 
the highest derivative term being </> itself, yielding (/> to be a function of u 
alone (note that x and t do not appear explicitly in any of the determining 
equations). Inserting this into the determining equations, one must have that 
is constant, a contradiction to our standing assumption in this subcase. 
Secondly, the process may terminate with an inconsistency, and thirdly, the 
process may terminate but with such a large expression that the result is useless. 

Case 3.2.4 ^ = u, /3 = 3, 7^ 0. Four determining equations were obtained, 
the first of which is = 0, so we substitute (j) = F{x,t)u + G{x,t) into 
the remaining three and require t) ^ 0. We find that there are no such 
solutions. 

Case 3.2.5^ = w, (5 = 3, 0„ = and not both n and a + 1 are zero. One 
determining equation was obtained which was a polynomial in u of degree two 
whose coefficients are functions oi x,t only, so the coefficients of powers of u 
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must be zero. These equations were easily simplified using the procedures in 
dif f grob2 to yield, 



K^Q, a = -l, (j) = 0, (3.9) 

K^O, a = -l, 0=—^, (3.10) 

t + Ci 

K arbitrary, a ^ —1, (f) — ci exp(C) + C2 exp(— ^), 

2Kt \ (3.11) 

1 + 



C — iy'a [ X — 



In (3.9) if we solve (1.1) and the invariant surface condition as a system of 
equati ons w e find that the only solution is u{x,t) = c, a constant. 

In (Islol) we can solve ( [l.l[) and the invariant surface condition to give the 
exact (canonical) solution 

u(x, t) — —2k + x/t, 

which cannot be realised by any of the previously found reductions, though 
it would not appear to be a particularly interesting solution. It is interesting 
to note that performing the KolRitt algorithm of diffgrob2 on the system 
comprising the original equation with the invariant surface condition led to 
a simple calculation for u. By contrast, the usual procedure of solving the 
invariant surface condition using the method of characteristics and inserting 
the result into the original equation to obtain the reduction was considerably 
more difficult due to the implicit nature of the reduction. 

In (3.11) we can again solve our problem to yield the exact (canonical) 
solution 

u{x, t) = ± (co + cie'- + C2e '') ' , <, = iv" 

I + a 

which is a special case of the travelling wave reduction 2.5. 

Case 3.2.6 — u, (3 — 3, (pu = 0, k = 0, a = —I. We are left simply with the 
determining equation (jj^x — = 0, which produces the following infinitesimal, 

(f, = g{t)e'= + h{t)e-'' , (3.12) 

where g and h are arbitrary functions. Hence we have to solve the invariant 
surface condition 

uux+ut^ git)e'' + h{t)e-''. (3.13) 

It is straightforward to show that every solution of this equation is also a solution 
of (P). 
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4 Nonclassical (r = 0) and Direct Methods 

In the canonical case of the nonclassical method when r = we set ^ = 1 
without loss of generality. We proceed by eliminating Ux,UxxiUxxx and Uxxt in 
( PI ) using the invariant surface condition ( ^.5|) which yields 

Ut - e(t)(j)uuUt - e0a;uWt " e^^U* - 4>xxU - (j)u4'xU - (j)'^(j)uuU - 2(t)(j)xuU 

(4-1) 

-(jxpiu - a4>u - (3(f)4>x - e0t0„ - (3(f)^4>u - e(f)xt - ^4>4'tu + 2k(/) = 0, 

which involves the infinitesimal </) that is to be determined. As in the t case 
we apply the classical Lie algorithm to this equation using the first prolongation 



pr'^-'^W and eliminate Ut using (4.1) 



The equivalent approach using the direct method of Clarkson and Kruskal 
[15] is to consider the ansatz u = U{x,t,w{t)) and require that the result be 
ordinary differential equation for w{t); see also [12,36]. It is straightforward to 
show that this yields the equivalent reductions. 

Case 4.1 e = 0. The nonclassical method generates a single equation of 25 
terms, without any singular solutions. Since this is difficult to solve explicitly, 
we seek polynomial solutions in u. 

Ansatz 1. (j) ^ F{x, t). In this case we obtain the following three exact solu- 



tions for (1.1) with e = 0: 

u{x, t) =^l2[x- (2k - (3ti)tf + /io, (4.2) 
where /i2 and /iq are arbitrary constants, provided that a ~ Q, 



u{x, t) = ^^^^ + ^.{x - 2Kt) + 6t^/'^, 
where 5 is an arbitrary constant, provided that a — Q and [3—^1, and 



(4.3) 



u[x,t)^ — , (4.4) 

at 

provided that a ^ 0. 

Ansatz 2. (p = F{x^t)v? + G{x,t)u + H{x,t). In this case we obtain the 



following three exact solutions for (1.1) with e ~ 0; 



u{x,t) ^ Ata,ii[^^/a{x - 2Kt)] , (4.5) 
where ^ is an arbitrary constant, provided that /3 = — 3, 

uix,t) = Aexp{n{x-2Kt)}, ^^2^-— (4.6) 

1 + p 

provided that /3 7^ — 1, and 

u{x,t) = Asech{^y/a{x - 2K.t)}, (4.7) 
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provided that /3 = — 3. 



Case 4.2 e = 1. In this case the nonclassical method generates a single 
equation of 150 terms, which has a singular solution if and only if 

(l>(l)u + 0^ - M - 2k//3 = 0, 

provided that a — f3 — 1 — 0. We again seek polynomial solutions of (p using 
one ansatz. 



Ansatz 1. 4> — F{x,t). In this case we obtain three following three exact 
solutions for ( |l.l| ) with e = 0: 

u{x, t) = ^2 [a; - (2k - (3^J)tf + /ii [x ~ (2k - I3^i)t] + /iq, (4.8) 
where /i2, /ii and /io are arbitrary constants, provided that a — 0, 

, , {x-2Ktf LL2{x-2Ktf fl + ^lA \, „ ^ ,1/2 

" i2t + ^ + + " ^ 

, M2(6 + 16^.1) , , ^ ^ 

H h2K + ^i +^2, 

(4.9) 

where /Li2, /xi and J are arbitrary constants, provided that a — Q and /3 = — 1, 
and 

X — 2nt 

u{x,t) = — , (4.10) 

at 

provided that a ^ 0. 



5 Discussion 

In this paper we have classified symmetry reductions of the nonlinear third order 
partial differential equation (1.1), which contains three special cases that have 
attracted considerable interest recently, using the classical Lie method and the 
nonclassical method due to Bluman and Cole [8] . The use of the MAPLE pack- 
age dif f grob2 was crucial in this classification procedure. In the classical case 
it identified the special cases of the parameters for which additional symmetries 
might occur whilst in the nonclassical case, the use of dif f grob2 rendered a 
daunting calculation tractable and thus solvable. 

In their recent paper, Gilson and Pickering [32] discuss the application of 
the Painleve tests for integrability due to Ablowitz, Ramani and Segur [3,4] 
and Weiss, Tabor and Carnevale [57] to equation (1.1). In particular, they 
investigate the integrability of the ordinary differential equations arising from 
the travelling-wave reductions 2.1 and 2.5 above. It would be interesting to in- 
vestigate the integrability of some of the ordinary differential equations arising 
from the other reductions derived in this paper using standard Painleve analy- 
sis, "weak Painleve analysis" [49,50] and "perturbative Painleve analysis" [20], 
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though we shall not pursue this further here. Marinakis and Bountis [41] have 
also applied Painleve analysis to the FFCH equation (1.4); an interesting aspect 
of their analysis is the use of a hodograph transformation. To conclude we re- 
mark that the RH equation (1.3) is a quasilinear partial differential equation of 
the form discussed by Clarkson, Fokas and Ablowitz [14]. It is routine to apply 
their algorithm, which involves a hodograph transformation, for applying the 
Painl eve PDE test to such quasilinear partial differential equations and show 
that (1^) does not satisfy the necessary conditions to be solvable by inverse 
scattering. 
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